Nuclear Excitation by a Zeptosecond Multi— MeV Laser Pulse 
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A zeptosecond multi-MeV laser pulse may either excite a "plasma" of strongly interacting nucleons 
or a collective mode. We derive the conditions on laser energy and photon number such that either of 
these scenarios is realized. We use the nuclear Giant Dipole Resonance as a representative example, 
and a random-matrix description of the fine-structure states and perturbation theory as tools. 

PACS numbers: 42.50Ct, 24.30Cz, 24.60Dr 
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Purpose. Qualitative Considerations. With the start 
of the construction of ELI (the "extreme light infrastruc- 
ture") Q or with extant ultra-intense laser facilites like 
NIF (if reconfigured as femtosecond pulse systems) @, HJ 
nuclear spectroscopy using intense high-energy laser 
beams with short pulses has become a realistic possibility. 
Indeed, it is envisaged to generate in the decade ahead 
pulsed laser light with photon energies of several MeV 
and pulse lengths of 10~ 19 seconds by coherent Thomson 
backscattering This will be possible provided that 

present intense experimental and theoretical efforts will 
validate the concept of an electron mirror Q . These very 
exciting developments call for a theoretical exploration of 
the expected nuclear excitation processes. In the frame- 
work of the nuclear shell model (a mean-field approach 
with a residual nucleon-nucleon interaction), two scenar- 
ioss come to mind, distinguished by time scales, (i) The 
time scale for the residual interaction is large compared 
to the time scale for laser excitation of individual nu- 
cleons. Then a single laser pulse containing N ^> I 
photons excites many nucleons more or less simultane- 
ously. The resulting "plasma" of excited and interacting 
nucleons (distantly similar to the initial stage of a pre- 
compound reaction) is instable. Nucleons excited above 
particle threshold with low angular momenta are emitted 
instantaneously. The remainder of the system is equili- 
brated by the residual interaction. Exciting questions 
are: What are the mean mass number, the mean exci- 
tation energy and the mean angular momentum of the 
resulting compound nucleus? How big are the spreads of 
these quantities? Presumably it will be possible to study 
compound nuclei at excitation energies and spin values 
not accessible so far. (ii) The time scale for the residual 
interaction is sufficiently short compared to the time be- 
tween two successive photon absorption processes. Then 
the nucleus relaxes after each photon absorption process. 
Single photon absorption leads to a collective mode (typ- 
ically the Giant Dipole Mode), and multiple photon ab- 
sorption within the same laser pulse may lead to the for- 
mation of higher harmonics of that mode. Thus, scenar- 
ios (i) and (ii) lead to extremely different forms of nuclear 
excitation. 

In this Letter we establish the time scales and the re- 
sulting conditions on the mean photon energy El and 



the number N of photons relevant for scenarios (i) and 
(ii). We do so by studying scenario (ii) in detail. We 
show that scenarios (i) and (ii) both occur for realistic 
choices of El and N. We also show that scenario (ii) is 
dominated by single-photon absorption. 

We focus attention on dipole absorption, the dominant 
photon absorption process in nuclei. The dipole mode 
1 10) is the normalized product of the dipole operator and 
the wave function |0) of the nuclear ground state. The 
dipole mode is not an eigenstate of the nuclear Hamilto- 
nian H nuc and is spread over the eigenstates \n) of H nuc 
with eigenvalues E^, fi = 1, Gross features versus ex- 
citation energy E of that spreading are measured by the 
strength function S(E) = \(10\n)\ 2 6{E - E^). The 
average (overbar) is taken over an energy interval large 
compared to the average nuclear level spacing d. In the 
simplest model adopted here, S(E) has Lorentzian shape 
and is characterized by two parameters 0] : The peak en- 
ergy i?di P ~ 80 A -1 / 3 MeV (where A is the nuclear mass), 
and the width T l sa 5 MeV (the "spreading width"). The 
resulting broad peak of S(E) is referred to as the Giant 
Dipole Resonance (GDR). By the uncertainty relation, 
the time for the dipole mode to spread over the eigen- 
states of -Hnuc (the "equilibration time") is r cq = ft/Fk 

The width Fdi p for gamma decay of the GDR to the 
nuclear ground state is estimated below and has a typ- 
ical value of 5 - 10 keV. With N coherent photons in 
the laser pulse, the characteristic time scale for pho- 
ton absorption is Tdi P = h/ (NT dip)- One would expect 
that excitation of the GDR (as opposed to scenario (i) 
considered above) dominates whenever Tdi P > r cq , i.e., 
whenever iVFdip < T^. That simple estimate is modi- 
fied by two factors, however, (i) For a short laser pulse 
with energy spread a (where we take a ~ 10 keV corre- 
sponding to a pulse length of « 10~ 19 s), the Lorentzian 
shape of the GDR produces for El < E d i P an addi- 
tional factor [T^/(El — Edi P )] 2 . (ii) The characteris- 
tic cubic dependence of Fdi p on energy yields an ad- 
ditional factor (El/ E^ip) 3 ■ In total the criterion for 
collective excitation of the GDR at energy El reads 
N < (E L /E di p) 3 [(EL-E di p) 2 /(T dip T±)}. With r dip = 10 
keV, r+ = 5 MeV, E dip = 14 MeV, E L = 7 MeV that 
yields N < 5 x 10 3 . That bound on N is significantly 
larger than the bound N < T^/T d i p sa 700 obtained from 
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the naive estimate and shows that even for an intense 
laser pulse, excitation of the collective GDR is a real- 
istic alternative in nuclei to multiple excitation of indi- 
vidual nucleons provided only that El is sufficiently far 
below -Edip- Thus, varying both N and El provides the 
exciting opportunity to investigate scenarios (i) and (ii) 
separately as well as the dynamical interplay between. 

We support these qualitative arguments by calculating 
the probabilities Pi for single-quantum dipole excitation 
and Pi for double-quantum dipole excitation as functions 
of N,El, and a. To calculate Pi we use the Brink- Axel 
hypothesis [1, Q . The hypothesis implies that single ex- 
citation of the dipole mode may be followed either by 
double excitation of that mode (i.e., formation of the sec- 
ond harmonic) or by dipole excitation of the configura- 
tions mixed with the single dipole mode. We account for 
both possibilities and show that for a -C the contri- 
bution from the Brink-Axel mechanism dominates and 
yields Pi = (l/2)P 1 2 . For values of N and of El such 
that Pi <C 1 that relation implies that single photon ab- 
sorption is the dominant process even if N 3> 1. Our 
result suggests that the probability for nuclear excita- 
tion by n-fold dipole absorption may be approximately 
given by P n « 2 _n P". That would imply that in the 
regime where our approximations apply (Pi < 1/2 or so) 
multiple collective nuclear excitation is unlikely. 

For the complex configurations that mix with the single 
or double dipole modes, we use a random-matrix model. 
Every such model is based upon the implicit assumption 
that the equilibration time (here r eq ) is short compared 
to the time scale of the physical process of interest (here 
Tdip)- Our use of random-matrix theory is justified if 
the above-mentioned conditions for collective excitation 
of the GDR are met. We also use perturbation theory 
to calculate Pi and Pi. That is justified if Pi and Pi 
are sufficiently small compared to unity. The resulting 
constraint is the same as for the use of random-matrix 
approach itself. 

Hamiltonian. We write the total Hamiltonian as 

H(t) = H nuc + H{t) (1) 

where H(t) stands for the time-dependent interaction 
with the laser pulse. In constructing H m , r we are guided 
by the following qualitative picture [l0(. In a closed- 
shell nucleus, the dipole mode 1 10) is a superposition of 
one-particle one-hole (lp lh) excitations. That mode 
is embedded in a sea of 2p 2h excitations |0fc) where 
k = 1 , . . . , K and K 3> 1 . (Here and in what follows 
the first label of the state vector counts the number of 
absorbed dipole quanta and the second enumerates the 
states). The mixing of both kinds of excitations causes 
the dipole mode to be distributed over the eigenstates of 
_ff nuc . The absorption of a second dipole quantum may 
either lead from the dipole mode 1 10) to the double dipole 
mode 1 20) (a 2p 2h state), or it may lead from one of the 
2p 2h states |0fc) to the dipole mode |lfc') of that same 



state (a 3p 3h state). The double dipole mode 1 20) is 
similarly embedded in a sea of 3p 3h states \0a) with 
a = 1, . . . ,L. All of the states \lk') are embedded in a 
sea of 4p 4h states \0p) where p — 1, . . . , M and M ^> K. 
The residual interaction of the nuclear shell model mixes 
these configurations, and both the double dipole mode 
and the states |l/c') are spread out over the eigenstates 
of -ff n uc- In modeling this qualitative picture we disregard 
the fact that single or double dipole excitation may pop- 
ulate states with different spin and isospin values. ]J nuc 
is accordingly schematically written in matrix form as 
follows. 
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Here Eq is the energy of the nuclear ground state, while 
Ei and E 2 are the mean excitation energies of the sin- 
gle and of the double dipole modes. For simplicity we 
use a harmonic-oscillator picture so that Ei — Ei = 
E\ — Eo = Edip. Moreover we put Eq = 0. The real 
matrix elements Vu mix the dipole mode with the 2p 
2h states 1 0/ ) . These are governed by the i^-dimensional 
Hamiltonian matrix . Similarly, the matrix elements 
Vip mix the double dipole mode with the 3p 3h states 
1 0/3) . These are governed by the L-dimensional Hamil- 
tonian matrix H^l. We write Hj^ — EiSki + and 

^a} = Ei5 a fj+H^p and assume that both H^]' and 
are random matrices, members of the Gaussian Orthog- 
onal Ensemble (GOE), with no correlations between the 
elements of Hj^' and of H^}. The spectra of EiSki+H^' 

(2) 

and of E 2 S a p + H^g both have the shape of a semicir- 
cle centered at Ei and E 2 , respectively. The last diag- 
onal block in Eq. @ describes similarly the mixing of 
the states \lk'} with the 4p 4h states |0p). We write 
Hk'l' = E 2 S k 'i' + T-Lk'i' and h a/3 = E 2 5 a p + h a p. We im- 
plement the Brink- Axel hypothesis by putting % = H^ 1 ' . 
Again, the M-dimensional matrix h pa is assumed to be a 
member of the GOE. We calculate the excitation proba- 
bilities Pi and Pi as ensemble averages for K, L, M —> oo. 
In that limit, the spreading widths of the single and 
double dipole mode and of each of the states |lfc') are 
given by the generic expression [llj = 2irv 2 p where v 2 
stands for the mean square of the relevant mixing matrix 
elements and p for the mean level density in the center 
of the semicircle. To avoid unneccessary complexity we 
assume that all spreading widths have the same value T^. 
That schematic picture can be refined if the need arises. 
We disregard the fact that the states excited by gamma 
absorption may decay by particle or by gamma emission. 
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That is justified because the time scales associated with 
such decay are orders of magnitude larger than both r eq 
and T dip . 

For the time-dependent interaction Hamiltonian H(t) 
we use a semiclassical description (justified for N 1) 
and write 



(3) 



H(t) = VNg(t)H di 



Here fldip is the time-independent electromagnetic in- 
teraction operator for a single-photon dipole transition. 
The factor y/~N accounts for the presence of TV ^ 1 pho- 
tons and the ensuing factor N in the transition rate. 
The dimensionless function g(t) describes the time de- 
pendence of the short laser pulse. We use the ansatz 



g(t) = exp[-crr/(2?r) - iu L t] 



(4) 



Fourier transformation of g(t) shows that the mean en- 
ergy of the laser pulse is — hu>L , the spread in energy 
has width a. Actually the interaction fldip depends on 
energy, too, via the wave number k. For a ss 10 keV we 
may put k sa fcj, where k^ — Ej,/(hc). 

In the scheme of Eq. ([2]) the non-zero matrix elements 
of the dipole operator are (10|fldi P |0), (20|fldi P |10), and 
( 1 fc'| fldip |0/c) . We use the Brink- Axel hypothesis to write 
(lk'\H dip \0k) = 5 kk , (lk\H dip \0k) . We assume that all 
non-zero matrix elements of the dipole operator have 
the same value written as (fldip). That corresponds to a 
harmonic-oscillator approximation. To estimate (fldip), 
we write the Hamiltonian H lni describing the interac- 
tion with the electromagnetic field in Coulomb gauge as 
flint = — (l/c)jA. Here j is the current and A the vector 
potential. In our time-dependent approach the latter has 
the form of a wave packet, 

A(x, Q,t) — a dui exp[— iuit] g(uj) exp[ifcr] x ■ (5) 



The unit vector x describes the polarization, 57 indicates 
the direction of the vector fc, and k = \f& and lo are 
related by k = ujje. The function g is the Fourier trans- 
form of g(t) in Eq. We determine the normalization 
constant a from the reqirement that the energy carried 
by A be equal to El. We use the dipole approxima- 
tion. That yields a 2 — (a El) /(it 1 / 2 he). Quantization 
of the electromagnetic field for individual quanta that 
have the form of the wave packet §E§ yields for the en- 
ergy density the expression n(E) = l/(47r 3 / 2 er). From 
Fermi's golden rule, the total width for dipole decay is 
r dip = 27m(£ L )|(fl d ip>| 2 . Thus, 



l(tfdip) 



2 7 r 1 /2r dip a 



(0) 



For the dipole width we use the Weisskopf estimate, 
r dip = li-(kR) 2 E L . With R = 3 x 1CT 13 cm and 
E L = 15 MeV that gives r dip w 10 keV, so that 



|(fldip)| ~ 10 keV, too. A somewhat larger value for 
Fdi p results when the Thomas-Reiche-Kuhn sum rule is 
taken into account. Here we are interested in order-of- 
magnitude estimates only, however. 

Perturbation Expansion. We solve the time-dependent 
Schrodinger equation in the interaction representation 
where the perturbation has the form 

H(t) = exp[iH nuc t/h] H(t) exp[-iH nuc t/H] . (7) 

We assume that at time t = — oo the nucleus is in the 
ground state |0). We determine perturbatively the prob- 
abilities Pi and P2 that at time t = +00 one or two dipole 
quanta have been absorbed. 

At t = +00, the probability amplitude for occupation 
of the state 1 10} reached after single-dipole absorption is 



61 = ^ 10 ' 



+ OO 



dt fl(t)|0) 



N 



ih 



(fldip) 



+ OO 



dt g(t)(10\exp[iH auc t/h]\10) (8) 



and analogously (with (10| in the last line replaced by 
(0fc|) for &ofc- The corresponding amplitudes for occupa- 
tion of the states |20), \0a) and |lfc'), \0p) reached after 
double-dipole absorption are denoted by 62, b a , b\k> and 
bo p . For example, we have 



bo, = ( E ) (OH 



N 
ih 



dh H(h) / dt 2 fl(t 2 )|0) 

J — OO 



fldi P r / d*i 9(h) / dt 2 g(t 2 ) 
^(0p|exp[-^fl nuc )^ 1 /n]|l^> ^ r 



xiOllexpiiiH^ih-^j/hWW) 



(9) 



The average probabilities for single and double dipole 
absorption are, thus, given by 



p? 



(10) 



The big angular brackets indicate the ensemble average. 
The first (last) two terms that contribute to P 2 are due to 
double excitation of the dipole mode and to the Brink 
Axel hypothesis, respectively. 

Averages. By way of example we perform the ensem- 
ble average for Pi and focus attention on the sum of 
the squares of the time-dependent matrix elements in 
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Eqs. Using completeness and a simple identity we 
obtain for these 



[10\exp[iH mu .(t 1 -t 2 )/K]\10)) = 
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We use Eq. © to write 
1 
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e - £ dip ± (*/2)r4. 



(12) 



Using Eq. (|4]) for g(i) and carrying out the time integrals 
(see Eqs. ((8])), we find that e is confined to an interval 
of size a around El- Since a -C T^, the argument of 
the expression in Eq. (fT2j) can be taken at e = El- The 
remaining integration can be done. With the help of 
Eq. © that yields 



Pi 



2^AT dip r+ 



(^L-Sdi P ) 2 + (i/4)(ri)2 



(13) 



The result (|T3|) is intuitively appealing and clearly dis- 
plays the suppression factors T^- /(El — E^p) 2 and 
(E'i/E'dip) 3 mentioned above that come into play for 
El < Edip- 

The calculation of P2 proceeds similarly but is more 
involved. We use operator identities such as 
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That leads to products of terms each containing H' 1 ' in 
the denominator. We neglect the correlations between 
eigenvalues of in different factors because such cor- 
relations extend over an energy range measured in units 
of the mean level spacing d while the range of the terms 
in Eq. (|14l) is given by ^> d. For the last two terms in 
the second of Eqs. (|10l) we obtain 
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with Pi given by Eq. (fT5|) . The calculation of the first 
two terms yields a contribution that in comparison to 



Eq. (|T5j) is small of order a/T^. Thus for all values of 
El the contribution to P2 from double excitation of the 
dipole mode is negligibly small in comparison with that 
from the Brink- Axel mechanism in Eq. (fl~5)) . As a result 
we find 



1 



-p? 



(16) 



The factor 1/2 in Eq. (fTrj)) is due to the time ordering 
in Eqs. Thus, we expect that for arbitrary positive 
integer n we have P n — 2~ n P™ . 

Conclusions. We have established the time scales and 
the resulting values for mean photon energy El and mean 
photon number N required for the realization of either 
of the two scenarios mentioned in the Introduction. This 
was done with the help of a random-matrix model for sce- 
nario (ii) for which we have calculated the probabilities 
Pi and P 2 for single and double nuclear dipole absorp- 
tion. Our assumptions and approximations require both 
Pi and P2 to be small compared to unity. In that case, 
scenario (ii) applies. Eq. (IT51) shows that in the tails of 
the GDR that condition is easily met even for an intense 
laser pulse. Ways of detecting such collective nuclear ex- 
citation experimentally are discussed in Ref . 12| . Double 
photon absorption is dominantly due to the Brink-Axel 
mechanism (as opposed to double excitation of the dipole 
mode). For Pi small compared to unity, single-photon 
absorption is the dominant mechanism while P2 <SC Pi ■ 

With increasing N, the time for dipole absorption Tdi p 
eventually becomes small compared to the nuclear equili- 
bration time T cq , and the competition between collective 
excitation and the formation of a strongly interacting nu- 
cleon "plasma" is decided in favor of the latter. Eq. (fl"3|) 
shows that in the center of the GDR (El ~ E d i p ), that 
will happen already for fairly small values of ./V w 10 or 
so. As N is increased, the process spreads to the tails 
of the GDR. It is a challenge to attain a theoretical un- 
derstanding of scenario (i), and of the interplay between 
both scenarios. 
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